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lo  Introduction 

It  has  been  observed  in  the  irrigation  wells  of  the  Hawaiian  Islands 
that  the  water-level  fluctuations  have  frequency  components  corresponding  to  those 
of  the  ocean  tides  QlI  .  This  phenomenon  was  analyzed  by  Carrier  and  Munk,  , 
assuming  the  observed  ground-water  fluctuations  to  represent  a  diffusive  trans¬ 
mission  of  the  tidal  disturbances  through  the  porous  volcanic  structure  of  the 
island.  The  purpose  of  the  investigation  was  to  use  the  results  in  estimating 
the  permeability  of  the  porous  medium. 

In  [2]  it  was  assumed  that  the  porous  medium  was  infinitely  deep.  In 
actual  fact,  however,  there  will  be  an  essentially  impenetrable  bounding  surface 
(see  Fig.  1).  This  paper  is  concerned  with  the  analysis  of  the  same  problem 
treated  in  [2]  ,  but  taking  account  of  the  bounding  bottom  surface.  Numerical 
computations  are  carried  out  for  several  values  of  the  dimensionless  depth.  Also 
the  limiting  case  of  shallow  water  theory  is  studied.  Using  the  results  of  the 
infinite  depth,  shallow  depth,  and  finite  depth  theory  it  is  possible  from  the 
graphs  given  in  Figs.  2  and  3  to  estimate  the  amplitude  and  phase  lag  in  the 
fluctuations  of  the  ground-water  as  a  function  of  the  distance  inland  for  various 
values  of  the  dimensionless  depth.  It  is  found  that  for  the  values  of  the 
physical  parameters  which  are  probably  of  most  concern  the  infinite  depth  theory 
gives  satisfactory  results  in  the  region  of  interest. 

2o  Formulation  of  the  Problem 

Although  the  formulation  and  the  first  part  of  the  analysis  of  this 
problem  follows  quite  closely  that  related  in  £2]  it  is  convenient,  for  the  sake 
of  completeness,  to  repeat  part  of  that  work  here.  As  the  water  in  the  ocean 


FIGURE  I 


(DOTTED  LINES  INDICATE  THE  EXPECTED  CORRECTION  TO  THE  COMPUTED  CURVES  WHEN 
ALL  MOOES  ARE  CONSIDERED) 
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bounded  by  0y<  and  OB  (see  Pig.  1)  rises  and  falls  about  its  mean  level  0Aa 

the  pressure  on  the  line  OB  varies.  Corresponding  to  this  periodic  change  in 

pressure  on  OB  we  can  expect  periodic  fluctuations  in  the  free  surface  of  the 
ground-water..  i,e«  Ox*. 

The  equations  governing  the  motion  of  the  fluid  in  the  porous  medium 
are  the  conservation  of  mass 

div  (pv)  -  -  ~  (ap),  (2.1) 

St 

and  Darcy's  law  which  replaces  the  conservation  of  momentum  law  (see  (3)  )a 

v  *  -  (k/p.)  grad  (p  -  po).  (2.2) 

Here  p  and  p  are  the  density  and  the  viscosity  of  the  fluids  &  and  k  are 
the  porosity  and  the  permeability  of  the  materials  v  is  the  velocity  with 
components  u  and  v  in  the  x5  and  y*  directions  p  is  the  gauge  pressures 

PQ  »  -  p0gy'  is  the  pressure  when  the  fluid  motion  is  zeros  PQ  is  the  mean 

densitys  and  subscript  notation  indicates  partial  differentiation.  The  simple 
compressibility  law  used  by  Carrier  and  Munk  is 

P«  =  Po'oC1  +  6  (P  -  P0)}  ^  (2*3) 

where  8  is  essentially  (p  c^)  with  c  the  speed  of  sound  in  the  fluid.  It 

o 

snould  be  pointed  out  that  Eq.  (2.2)  says  that  the  pressure  gradient  is  . 

proportional  to  a  velocity  rather  than  an  acceleration  as  in  the  Navier-Stokes 

equation.  As  a  result  we  will  obtain  finally  an  equation  of  the  diffusion  type 

rather  than  a  wave  equations  hence  the  free  surface  amplitude  will  decay  in  x*. 

The  boundary  conditions  expressed  in  terms  of  the  pressure  p  are 

p  *»  0,  on  the  free  surface  (2«lia) 

P  ■  -  pQgy'og  ♦  on  OB  (2.iib) 

9(p  -  P  )/dy'  “0  on  y*  *  <=  H,  (2.Uc) 

o 

where  H  is  the  depth  of  the  porous  medium.  The  last  boundary  condition  results 


from  the  requirement  that  the  normal  component  pf  velocity  be  zero  on  the 
impenetrable  bottom.  The  pressure  is  of  c purse  directly  proportipnal  to  the 
tidal -wave  amplitude. 

If  we  let  q  -  p  -  PqS  and  combine  Eqs.  (2.1),  (2.2)  and  (2.3)  we 

obtain 

at  Q 

qt,  (2.3) 

where  A  is  the  Laplacian  operator.  If  we  denote  by  r)(x*,t)  the  y5 
coordinate  of  the  free  surface,  Eq.  (2.4a)  implies  q(x,Jrj,t)  =  pQgr),  but  on  the 
free  surface  *  v/<9  ,  hence  using  Eq.  (2.2)  our  boundary  condition  (2.4a)  may 
be  expressed  as 

qt  +  qy.,  *»  0  on  y»  **  tj,  x  >  0.  (2.6a) 

^  o 

Actually  as  in  the  usual  linear  theory  of  water  waves  this  boundary  condition  is 
to  be  applied  on  y  =  0.  The  boundary  conditions  (2.4b)  and  (2.4c)  may  be 
written  as 

q  (on  OB)  =  q^eicot,  (2.6b) 

■  0,  y»  »  -  H.  (2.6c) 

We  shall  only  solve  this  problem- in  the  case  that  the  line  OB  occupies  the 
half°line  y  ®  0,  x «  0.  That  is  we  take  £  =  0°  (see  Fig.  l).  Actually  this 
is  fairly  realistic  since  %  is  probably  of  the  order  of  3°  or  so. 

Finally  if  we  introduce  the  following  dimensionless  variables 
trs  cot,  x  =  x!/Lg  y  =  y'/Lj  h  =  HZ!!  , 

L  s  (p0gk)/(p$oco),  e  »  (Pog2ks)/(ti0oco),  q  =  q1^(x,y)e1“t  , 

we  obtain 

A ((>  ~  Uf  ■  0,  .  (2.8) 


with  the  boundary  conditions 
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(2.9a) 

(2.9b) 

(2.9c) 


Cf>, '  +  1(J>  a  0  ,  y  “  0  ,  X  >  0 

V 

<£>  =  0  ,  y  =  -h  s  -  o©  <  x  <  00  . 

7 

-1 

The  free  surface  T)(x,t)  is  (pQg)  q(x,o,t)  but  from  Eqs.  (2.7) 
q  »  q  t/>  (xpy)  exp  (i©t)  so 

r](x,t)  *»  (x,o)eicot  x  >  0.  (2.10) 

P0S 

So  the  problem  of  determining  the  free  surface  is  exactly  that  of  determining 
p  (x,o).  The  combination  of  parameters  q^/pQg  is  the  maximum  height  of  the 
tidal -wave  measured  from  y 5  =0. 

Before  proceeding  to  a  solution  of  the  problem  defined  by  Eqs.  (2.8) 

and  (2.9)  it  is  perhaps  worth-while  to  mention  briefly  the  size  of  the 

-2  2 

parameters  which  appear  in  this  problem.  "We  have  p/p  =  0(10  cm  /sec), 

v  0 

“O  2  2 

k  ~  0(5x10  cm  ),  &  -■  0(,20),  g  «  980  cm/sec  ,  and  co  for  a  twenty  four  hour 
tide  is  2n/2l*  hours,  hence  L  =  0(1000  ft).  Since  c  is  0(5000  ft/sec) 
for  water,  s  is  a  very  small  number,  0(10"^).  Finally,  a  reasonable  value  for 
the  depth  of  the  ocean  is  about  three  miles  so  h  may  be  as  large  as  15. 

•; .  Shallow  Water  Theory 

Before  considering  the  general  problem  given  by  Eqs.  (2,8)  and  (2.9) 
let  us  look  at  the  limiting  case  in  which  the  depth  H  is  small  enough  that  we 
can  neglect  variations  in  the  y!  direction,  and  also  set  v  s  o.  Then 
u(x’ ,t)  represents  an  averaged  velocity  across  the  section  -H<y»<0.  If  we 
assume  incompressibility,  i.e.  5  =  0,  the  conservation  of  mass  equation 
appropriate  to  this  situation  is 

HP0V  =  "  Po^t  • 


Darcy ” s  equation  reduces  to 


(3.1) 


■  -s- 


u  »  -  (k/n)q  .  (3.2) 

x* 

Since  there  is  no  variation  in  the  y  direction  onr  condition  that  q  ■  pQgr] 
on  the  free  surface  must  hold  throughout  the  strip  — H<y'<0,  -x  >  0.  Using  this 
and  Eqs.  (3.1)  and  (3.2)  we  obtain 


^•*•-555*  ■0*  x‘  * °- 

The  condition  that  q  *  q  exp  (loot)  for  x  0  is  now  applied  at  x  ■  0$ 
consequently  we  set  q  »  q  $»(x‘)  exp  (ioat).  Then  Eq«  (3.3)  becomes 

i(HL)”1  0,  x«*0. 

An  appropriate  solution  of  Eq.  (3«U)  satisfying  a  finiteness  condition  at 
infinity  is 

(xs  jt)  *  e 

Hence 


(3.3) 


(3.U) 


x* 


ri(x«st)  =  (q-j/p0g)  e  ^2HL^ 

x 

■  (d/p^)  e  ^ 


1/2 


~(i/HL)l/2x* 

if  cab  -  -2Ei_  1/2\ 

•  l  (2HL)1/2J 

*  > 

(3.3a) 

1  . 

(3.3b) 

Actually,  in  order  for  this  theory  to  be  valid  not  only  must  the  wave 
length  of  the  disturbance  be  large  compared  to  H  as  in  the  usual  shallow  water 
theory  but  also  H  must  be  small  compared  to  the  other  natural  length  scale,  L, 
which  appears  in  the  problem,  i.e.  h  must  be  small.  This  can  be  seen  by  an 
examination  of  the  behaviour  of  the  solution  of  the  general  problem.  This  is 
done  in  Sec.  3  where  it  is  found  that  for  h  $  l/U  we  can  except  the  shallow 
water  theory  to  be  quite  accurate.  The  amplitude  and  phase  lag.  of  p  gT)(x,t)/q 
are  plotted  in  Figs.  2  and  3  as  a  function  of  x  for  h  =  l/U. 

Uo  Solution  of  the  Problem 

To  solve  the  problem  defined  by  Eqs,  (2.8)  and  (2.9)  we  shall  use  the 
method  of  Fourier  transforms  and,  the  Wiener-Hopf  technique.  Let 


where 
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clear  from  Eq.  (H*W  that  our  problem  is  now  that  of  determining  g^(x)  and 

hence  g  (  2).  To  determine  g  (?)  using  Eq.  (li.8)  we  shall  use  the  Wiener-Hopf 
2  2 

technique.  This  technique  has  been  used  to  treat  similar  problems  (see  for 
instance  Q>,  U,  5])  5  consequently  the  analysis  will  only  be  briefly  outlined 
here.  First  g_^(  S  )  is  analytic  in  the  upper  half  plane,  (UHP),  Im(  ?  )  ?»  -  a; 
gg(  %)  is  analytic  in  the  LHP  and  f(  ?  )  is  analytic  in  the  UHP.  The  function 
G( ?  )  is  analytic  and  non-vanishing  in  a  strip  containing  the  real  axis.  This 
will  be  seen  clearly  at  the  end  of  this  section  where  G(  E  )  is  represented  as 
the  quotient  of  two  infinite  products.  It  might  be  noted  that  though 
C  ~  (§  2  +  is)^2  is  a  multivalued  function,  G(S  )  as  defined  by  Eq.  (ii.9) 
consists  only  of  even  terms  and  hence  does  not  have  any  branch  points.  Assuming 
for  the  moment  that  we  can  write  G(  ?  )  as  G  (?  )/G+(  E  )  where  G  (  %)  is 
analytic  and  npn-vanishing  in  the  LHP,  and  G  (  ?  )  is  analytic  and  non-Yanishing 
in  the  UHP  we  can  rewrite  Eq.  (U.8)  as 


f(«  )  Q+(  «)  -  Gj-ia)  )  =  j\U)  -  Oj-ia)}  g^g  ) 


♦  0J5)  g2(«). 


(lulQ) 


The  left  hand  side  of  this  equation  is  analytic  in  the  UHP,  the  right  hand  side 
is  analytic  in  the  LHP  and  they  agree  in  a  common  strip  of  analyticity.  Hence 
Eq.  (Iwl0)  defines  an  entire  function  E( g ).  We  shall  show  shortly  that 

G  (£)  -  0(%V2)  as  |?  |  — *oo  ,  Im(g  )  <  0  and  G  (1*  )  =  0(£  -1/2)  as 

€=»  + 

|? |  — ^ oo  s  Im(g  )  >  0.  Using  this  and  investigating  order  conditions  at  infinity 
we  can  show  that  E(g)  **  0,  consequently 


g,(.S)  -  /B-(-ia)  -l\  5,(4). 
2  •  8.(4)- 


(lull) 


It  is  now  necessary'-  to  determine  G  ( E  )  and  G  ( %  )„  The  splitting 
of  G(  2  )  is  done  in  a  manner  exactly  analogous  to  that  used  by  Heins  in  [kj 
and  [$]  «  Using  the  infinite  product  representation  of  cosh  z,  (pee  U> ])  we 
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have 


where 


oo  \ 

M(2  )  *  cosh  Ch  =»  Tf  fi  +  Up2h2/(2n  +  l)2n2/=  m(fg  )  m(-4),  (U.12) 

n^o  v  ' 

"»  1/2 


-“>  •  m  •  agy 


n-o 


k  i2hg 
(2n+l)n 


-  igk  % 

e  (2n*l)n 


(U.13) 


We  have  inserted  the  exponentials  to  insure  absolute  convergence  of  the  infinite 
products  defining  m(  S  )  and  m(-  2  )  in  the  UHP  and  LHP  respectively.  If  we 
■write  M(  $  )  as  M  (g  )/M+ (4  )  .  it  is .  clear-  that  M(g)=m(g)  has  no  zeros 
or  poles  in  the  LHP.  Similarly  l/M+(  g  )  ■  m(-  4)  has  no  zeros  or  poles  in  the 
UHP. 

The  function  L( %  )  ■  C  sinh  Ch  +  i  cosh  Ch  has  zeros  at 
Ch  *»  *  ip^j  n-0,1,.2,  where  the  {3^  are  complex  numbers  lying  in  the  first 

quadrant.  For  n  large  they  may  be  determined  by  the  asymptotic  relation 
=■>  nn  +  ih/nn  +  0(  QnnJ’2).  We  may  write  L(  B  )  as 

oo 

.71 


where 


L(«) 


1(g) 


1  ♦ 


C2h2 


n=o 


vl/2 


il(4)  l(-£). 


(ii.ll;) 


1  + 


ich 

P? 


Again  we  have  inserted  the  exponentials  in  order  to  insure  absolute  convergence 
in  the  appropriate  half  planes.  If  we  write  L(  £  )  as  L  ( %  )/L  (  g  )  and  take 
Lj  g)  =  1(  4)s  l/L+(  g)  =  il(-g)  it  is  clear  that  L_(  4  )  is  free  of  zeros 
and  poles  in  the  LHP  and  L  (  g )  is  free  of  zeros  and  poles  in  the  UHP. 


Consequently  we  have 

X(*). 


0(8) 


V.s) 


X(g) 


n.(s ) 
AT(g) 


e  1(S) 

m(  £)  r 


L+(%)  _  e^m(-g) 


Mjg) 


il(-  £) 


(l;,l6a) 

(U.l6b) 


i 


We  shall  choose  the  factor  exp  (*<&)]  introduced  in  Eqs.  (i|.l6a)  and  (it.l6b) 
in  such  a  manner  that  G  (£;  )  and  G  ( <£  )  have  algebraic  behaviour  as 

M  + 

|gj  — I,  qo  in  the  LHP  and  OHP  respectively. 

To  investigate  the  behaviour  of  G_(  £  )  for  Im(  §  )  <  0,  |£|— *■ 
we  first  note  that  the  terms  involving  e  may  be  neglected  against  unity  for 
jgf  ~+oo  0  Since  |3  njt  as  n  — ♦  we  have  that  G_( S,  )  is  of  the  order 

e^(S)(l  *  iSh/8  )  ft  (1  ♦  */n)e~W/n  «  e*(^(l  +  IS  h/p  )(lT  (1  +  w/n)e“W/n 


7T/l  +  2w/(2n+l)l 

nssn  >*  ■/ 


-2w/(2n+l] 


7 T  (i  +  2w/n)  e 

n=l 


.  ln=l 

--iv7n 


where  w  =  i£h/n.  Now  using  the  relation  that 

yw  o°  -w/n 

1/P  (w)  =  w  e  77"  (1  +  w/n)e  s 
n=l 

and  Stirling's  asymptotic  formula  for  the  gamma  function,  (see  [6])  we  obtain 

Go(  2  )  ■  0  jV^exp  [X  (  S)  +  w  In  > 

for  Im(g)  <  0,  l&l'-+co  t  So  choosing  ^  (S)  ■  -  w  In  «  -  (igh  l11  U)/n  we 

have  that  Gjz)  -  0(2 1/2)  as  ,  Im(2  )  <  0.  A  similar  argument  will  show 

that  G+(2)  =  0(w“1/2)  for  |5|— 5  Im(2  )  S'  0.  With  these  order  relations  it  is 
not  difficult  to  show  that  E(  g- )  is  zero  as  mentioned  earlier^  and  hence  we 
obtain  Eq.  (ii.ll)  for  g^g)  where  G  (g)  is  defined  by  Eq.  (U.l6a).  In 
particular  it  can  be  seen  from  Eq.  (i;.l6a)  and  the  definitions  of  1(g)  and 
m(g)  that  G  (0)  =  1. 


Using  the  usual  inversion  formula  we  have  that 

OO  . 

g-(x)  »  ~  ll  d  £  s 

2  2"J  jcL(ft)  ]TS3?) 

-o°  OO  >»  ^ 

a  1  O-(-ia)  cosh  Ch  -  G+(g)G(g)  d  £ 


G+(2)G(^; 


a-i  g 


(U.17a) 


(lulTb) 


-  OO 
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In  the  limit  as  a  — ►  0  it  is  clear  from  Eq.  (lj.,17a),  that  g^(4  )  is  not 
singular  at  the  origin 5  hence  we  may  actually  take  the  real  axis  as  our  path  of 
integration  in  evaluating  g2(x)«  Of  course  we  shall  actually  close  our  path  of 
integration  in  the  UHP  when  x  >  0  and  in  the  LHP  when  x  <  0.  For  the  case 
x  >  0  it  is  convenient  to  use  Eq0  (Ij.lTb)  for  evaluating  g^(x).  Since 
G  (4  )  is  non-vanishing  in  the  UHP,  g  (x)  will  be  simply  the  sum  of  the 
residues  at  the  poles  of  the  integrand  which  occur  at  Ch  «*  ip^,  n  =  0,  l,"00*. 
Carrying  out  this  straight  forward  computation  we  obtain  in  the  limit  a  *  0 


2  2  1 /* 

where  =>  (p£/h  +  is)  «  In  the  particular  case  that  the  fluid  is  in¬ 
compressible,  i.e.  e  *  0,  Pn/an  ■  h  and  Eq.  (iul8)  becomes 


5.  Numerical  Computations  and  Discussion 

In  this  section  we  shall  only  be  concerned  with  the  case  when  the  fluid 

may  be  considered  as  incompressible,  then  g^(x)  is  given  by  Eq.  (H.19).  In  [2] 

a  few  values  of  g^(x)  were  computed  for  s  =  »01  and  compared  to  the  6  =  0 

case§  the  amplitude  and  phase  lag  in  the  ground-water  fluctuation  for  »  =  .01 

were  slightly  lower  than  for  6  =  0.  Since  e  is  however  0(10“^)  we  should 

expect  very  little  error  in  actually  setting  6=0. 

First  let  us  determine  when  the  shallow  water  theory  solution  given  by 

Eqs.  (3o5)  may  be  expected  to  be  valid.  In  order  for  g^(x)  as  given  by 

Eq.  (I4.0 1?)  to  agree  with  the  shallow  water  solution  it  is  necessary  that 
l/2 

aQ—  (i/h)  -as  h  — ♦  0,  and  also  that  all  the  coefficients  of  the  higher 


order  terms'*'  must  approach  zero.  It  can  be  shown  with  little  difficulty  from  an 

Investigation  of  the  transcendental  equation  C  sinh  Ch  +  i  cosh  Ch  that  for 

small  h,  (ih)'^  and  hence  a^*-  Also  upon  noting  that  for  small 

h,  nn  for  n  ^  1,  it  can  be  seen  from  Eq.  (U.19)  with  the  aid  of  the 

representation  of  G  (ict  )  given  in  Eq.  (5.1)  that  all  the  coefficients  of  the 

41  n 

higher  order  terms  do  approach  zero  as  h-p*  0.  Hence  gg(x)  as  given  by 
Eq.  (I4.0I.9)  does  approach  the  shallow  water  solution  as  h  0.  To  determine 
quantitatively  when  Eq.  (3.5b)  is  valid  we  have  computed  |3q  and  aQ  as  a 
function  of  h.  Also  the  ratios  of  the  wave  length  predicted  by  the  shallow 
water  theory,  X  =  2n/(2EL)^^  s  to  H  and  that  of  the  fundamental  mode, 

Xq  «  2n/lm(a0),  to  H  have  been  computed.  These  results  are  given  in  Table  1 
and  graphically  in  Fig.  Ij..  It  appears  from  Fig.  it  that  we  may  expect  the 
shallow  water  theory  to  be  accurate  over  the  entire  range  of  x  for  h  <  l/lt. 

In  order  to  compute  g  (x)  for  various  values  of  h  it  is  necessary 
to  cast  G+(ian)  into  a  form  more  suitable  for  numerical  analysis  than  that 
given  by  Eq.  (U.l6a).  This  can  be  done  in  a  straight  forward  manner  by  using 
the  infinite  product  representation  of  the  gamma  function.  We  obtain,  when 
e  *»  0,  that 

G+(i«n)  •■=  ?n  e*P  W») 

2ni(l  + 

In  any  actual  numerical  computation  the  infinite  products  in  Eq.  (5.1)  is,  of 

course,  to  be  replaced  by  a  finite  numbers  of  terms  (recall  that  8  — *  mn  as 

m 

m  cc).  The  number  of  terms  that  is  required  to  give  an  accurate  answer  is  of 
course  dependent  upon  h  and  Pn. 

\fe  shall  refer  rto  the  term  exp(~aQx5  in  Eq.  (U.lp),  w}iich  is  the  dominating 
term  as  x  =-►  00  ,  as  the  fundamental  term. 


(t<Pi/*)}2  fff  (1  *  MU 

r(2pn/n)  (rn*l  (1  +  Pn/mn) 


-1 


(5.1) 


r 


i 

%  1 


FIG.  4  DIMENSIONLESS  WAVE  NUMBER  vs.  h 

X  IS  THE  WAVE  LENGTH  PREDICTED  BY  SHALLOW  WATER  THEORY. 
X0  IS  THE  WAVE  LENGTH  OF  THE  FUNDAMENTAL  MODE  PREDICTED 
BY  FINITE  BOTTOM  THEORY. 
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In  this  paper  g^(x)  was  computed  for  h  ®  2  and  h  =  15.  Values  of 
pn  for  h  a  2  and  h  =  15  are  given  in  Table  2.  Let  us  first  consider  the  case 
h  =  2.  Then  the  real  parts  of  Cq  and  are  given  by  .59  and  1.65  respectively 
hence  we  may  expect  the  fundamental  term  to  give  an  accurate  result  for  x  2. 
Carrying  out  the  necessary  computations  we  obtain 


£°£r)(x,t)  -  .61  e“’59x  ei(®t-.885~.292x)  +  0(e”fllX)|  h  -  2 
<31 


(5.2) 


In  Figs.  3  and  U  the  amplitude  and  phase  lag  in  the  ground-water  fluctuation  have 
been  plotted.  The  extrapolation  of  the  results  to  x  =  0  are  indicated  by  dotted 
lines.  Actually  it  is  not  difficult  to  obtain  another  term  in  the  series,  but 
unless  particular  quantitative  information  is  desired  for  this  value  of  h  it 
hardly  seems  necessary  to  do  that.  It  might  be  mentioned  that  three  terms  were 
more  than  sufficient  in  evaluating  the  infinite  product  in  Eq.  (5.1). 

In  the  case  that  h'«  15,  the  fundamental  term  can  only  be  expected  to 
be  accurate  for  x  ^  7.  We  obtain 

pog  (  j_  \  n<7i  -.lOlpc  i(cot-1.68-.007x)  .  -a-jx 

Ji-  T)(x,t)  «  .G7i+  e  e  +  0(e  )$  h  -  15.  (5.3) 


In  order  to  obtain  results  valid  for  x  =  1  or  2  when  h  =  l5  would  probably 
require  the  computation  of  three  or  four  terms  of  the  series.  However  in  view 
of  the  results  for  h  =  2  and  these  results  for  x  7  it  is  clear  that  the 
amplitude  curve  for  h  =  15  will  lie  almost  exactly  on  the  curve  given  by  the 

O 

infinite  depth  theory^ (see  Fig.  2). 

It  is  interesting  to  note  that  the  shallow  water  theory  and  finite 
depth  theory  predict  an  exponential  decay  in  x  for  the  amplitude  of  the  ground- 
water  fluctuation!  this  is  in  contrast  to  the  algebraic  decay,  (like  x”^), 


2 

The  amplitude  and  phase  lpg  curves  for  h  =  oo  have  been  taken  from  the  results 
given  in  [2]. 
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predicted  by  the  infinite  depth  -theory.  Also  the  phase  lag  predicted  by  the 
shallow  water  theory  and  finite  depth  theory  continues  to  increase  with  x,  while 
that  predicted  by  the  infinite  depth  theory  approaches  n/2  with  increasing  x. 
(This  is  clearly  illustrated  in  Fig.  3).  That  gg(x)  as  given  by  Eq.  (h.19) 
approaches  the  infinite  depth  result  given  in  [2]  as  h-»<»  ,  cannot  be  seen 
easily  from  (h«19)o  However  an  examination  of  G(  4  )  as  given  by  Eq.  (In  9) 
shows  that  as  Ch  -*» o  9  G(£)  -  i  +  (&£  +  which  we  might  denote  by 

)  „  This  function  is  the  one  that  occurs  in  [2] .  It  is  interesting  to 
note  that  )  has  singularities  of  the  branch  point  type,  and  in  the  limit 

as  ®  ->  0  these  singularities  will  occur  at  the  origin.  This  explains  the 
algebraic  bshaviour  of  n(x,t)  for  h  =  at)  .  In  contrast  for  any  finite  value  of 
h  the  strip  of  analytic ity  of  G(&  )  is  finite  even  when  ®  -*  0,  and  its 
singularities  are  poles  rather  than  branch  points ;  hence  the  exponential  sort  of 
behaviour  for  r)(x,t)  for  finite  h. 

A  plausible  physical  explanation  for  the  fact  that  the  amplitude  curves 
fop  the  ground-water  fluctuation  lie  continuously  below  one  another  as  h 
depreases  (see  Fig.  2)  is  the  following.  Imagine  that  our  porous  medium  and 
fluid  occupy  the  strip  -H<y<0,  -  °°  <  x  <00  .  Suppose  that  we  apply  a 
uniform  pressure  on  the  half  line  y  »  0,  =  00  <  x  <  0;  then  fluid  in  the  left 
half  strip  will  be  forced  through  the  gap  -  h  <  y  <  0  and  the  free  surface 
given  originally  by  y  =  0,  x^O  will  rise.  The  amount  of  fluid  that  can  be 
forced  through  this  gap,  and  hence  the  effect  that  the  pressure  variation  can 
have  on  the  free  surface,  is  proportional  to  the  gap  distance,  h.  So  with 
decreasing  h  the  amplitude  of  the  free  surface  fluctuation  is  lower  and  dies 
out  more  quickly. 


The  author  would  like  to  express  his  appreciation  to  Professor 
G.  F„  Carrier  for  suggesting  this,  problem  and  for  his  continued  interest 
during  several  stimulating  discussions.  He  is  also  indebted  to  Professor 
R.  E.  Kronauer  for  his  helpful  comments. 


h 

Po. 

ao(e*0) 

X/H 

VH 

<=>*►  0 
.25 

- 

,36%  *  i  .3382 

1.1*70|*  +  i  1-3528 

17.76 

18.57 

-50 

.5376  *  i  .1*51*8 

1.0752  *.  i  .9096 

12.57 

13.81 

1-00 

.8001*  +  i  -5702 

.800}*  +  i  -5702 

8,88 

11.09 

2-00 

1.1828  +  i  -5832 

.5911*  +  i  .2916 

6.28 

10.77 

3-00 

1.3739  +  i  .1*775 

.1*579  +  i  -1592 

1*.85 

13.16 

5-00 

1-5033  +  i  .3090 

-3007  +  i  .0618 

3.97 

20.03 

10-00 

1.55U7  +  i  -1569 

.1555  *  i  .0157 

2.81 

1*0.01* 

15.00 , 

1.5638  +  i  .101*6 

.101*5  +  i  .0070 

2.29 

59.83 

TABLE  1. 
l/2 

h  is  the  dimensionless  depths  X  =  2s(2HL)  '  is  the  wave  length  for  shallow 
water  theory,  Xq  »  2nL/lm(ao3  is  the  wave  length  of  the  fundamental  mode  for 
finite  bottom  theory* 


O  H  t\i  (r\_=f'Lr\\o 


n 

Pn(h»2) 

Pn(h  »  15) 

•  0 

1 

2 

3 

1* 

3 

6 

1.1828  +  i  .5832 
3.3106  +  i  .61*99 
6.3011*  +  i  .3277 
9.1*21*8  +  i  .2138 

1.5638  +  i  .101*6 

1**6886  +  1  .3231* 

7.8016  +  ,571*9 

10,8705  +  i  .9070 

13,7139  +  i  1.3629 
16,1332  +  i  1.1*236 
18.961*1*  +  i  1.0572 

TABLE  2. 
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